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MOLECULAR VIBRATION SPECTRA FROM FIELD EMISSION 
ENERGY DISTRIBUTIONS 
by Dennis J. Flood 
Lewis Research Center 

SUMMARY 

The effect on a field emission current of foreign molecules located at the cathode 
surface is presented. The electron-molecule interaction is represented by the potential 
energy of a charged particle in the presence of an electric dipole field. Only the case of 
a quadratic dispersion law for the electrons in the metal is considered. Results have 
been obtained which are valid for the low temperature range where the tail of the Fermi 
distribution function may be neglected. The calculation indicates that the emission cur- 
rent should change by about 1 percent for less than monolayer coverages of the cathode. 
The total energy distribution of the emitted electrons P T (e), and in particular the first 
derivative of P T (e), should give information on the vibrational spectra of the molecular 
species involved. 


INTRODUCTION 

It has recently been pointed out (refs. 1 and 2) that vibrational excitations of foreign 
molecules located in the oxide layer of thin film metal-oxide- metal (MOM) sandwiches 
produce considerable structure in the current- voltage characteristics of the tunneling 
current. The effect shows up as a series of peaks in d I/dV . The peaks appear be- 
cause a fraction of tunneling electrons interact inelastically with the molecules, exciting 
them into one or another of their upper vibrational states. A model calculation (ref. 3), 
in which the electron-molecule interaction is taken into account by including the molecu- 
lar dipole field in the potential describing the barrier region, is in good agreement with 
experimental observation. This report presents a calculation based on the same kind of 
model for the case of field emission into a vacuum from a metal with a molecular layer 
on its surface. 



DESCRIPTION OF MODEL 


The field emission process from clean metallic surfaces can be treated by consider- 
ing a one-dimensional potential which is assumed to represent the effects of the metal 
and applied electric field on the conduction electrons. In the free-electron approxima- 
tion the effective potential energy of the electrons is written as 


inside the metal, and 


V = W„ - eFx - 
o a 


4x 


( 2 ) 


outside. Figure 1 is an illustration of the form of the effective potential energy. The 



Figure 1. -Effective potential energy of an electron near a metallic surface in vacuum, 
under the influence of a strong electric field. 


origin of the three terms in equation (2) is as follows: 

(1) W is the height of the potential energy well representing the metal. Energies 
will be measured from the bottom of the conduction band. 

(2) F is the applied electric field and is assumed to be constant as far as surface 
barrier effects are concerned. The origin of coordinates is chosen to be at the metal- 
vacuum interface. Positive values of x are measured to the right of, and perpendicu- 
lar to, the surface of the metal. Under these conditions an electron of charge -e out- 
side the metal will have a potential energy of -eFx. 

(3) The last term arises from the force exerted on the electron by the image charge 
induced on the surface of the metal as the electron attempts to escape. With the zero 


2 



of energy chosen as indicated, in the absence of an applied electric field an electron at 
infinity will have an energy W„ . 

Figure 1 also shows the classically forbidden region where the potential energy is 
greater than the kinetic energy of the electrons. The points x & and x^ are the classi- 
cal turning points of the motion. The presence of molecules in the barrier region may be 
taken into account by including in equation (2) a term arising from the Coulomb potential 
of the molecules and their images. Approximating this potential by a dipole potential, 
and neglecting the transverse components (i. e. , those parallel to the metal surface), the 
electron-molecule interaction may be represented by 


V t (x) = 


-2eP x x 



(3) 


The effect of this term will be to lower slightly the barrier height in the region where the 
molecule is located. The quantity P x is the component of the dipole moment normal to 
the surface, and r^ is the distance in the plane of the barrier of the electron from the 
molecule. The other components of the dipole field of the molecule and its image effec- 
tively cancel each other at distances greater than their separation. This is shown in the 
appendix. 


Electron Tunneling 


The probability per unit time that an electron in a state 
will make a transition to a state b on the other side is (ref. 


P ab- 


2ir 


M 


ab 


(W 1 - V 


a on one side of a barrier 
4) 


(4) 


where is the density of states at b, and f & and f^ are the probabilities of occupa- 
tion of the states a and b, respectively. The quantity is the tunneling matrix 

element for a transition from a to b. Bardeen (ref. 5) has shown that the matrix ele- 
ment may be written as 


M ab = ~ iKJ ib 


(5) 
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where is the matrix element of the x component of the current density operator 
J , evaluated in the barrier region: 


J ib = <^al J x|f'b> »> 

and and xjy^ are the wave functions which describe the states a and b. van- 

ishes unless the transverse wave number k|. is conserved in the tunneling process. If 
state a is considered to be on the left of a barrier, the current flowing from left to 
right is obtained by summing P ^ over all states a of fixed k t , summing over k t , 
multiplying by the electronic charge e, and multiplying by 2 to account for spin. Simi- 
larly, 

p ba = f KalVb^-V (7) 

is the probability per unit time for a transition from the state b to the state a. The 
net current density flowing from left to right across the barrier is thus given by 


J = 2e 


If 


( P ab p a " P ba p b^ de 


4?re 

fi 



(i - f b >p 


a 


M ba| 2p a f b (1 - f a )p b] d€ 


(8) 


The density of states factors allows the sum over the initial states to be replaced by an 
integral. But since J is equal to its complex conjugate J* 

A X 

M ba = ' 1R <*bl J xk>=- M ab «*> 



(9b) 


and equation (8) reduces to 



r 


•> = /j M ab| 2 <W f a " V d£ < 10 > 

where e is the total energy for a fixed k^. For the case of field emission from a metal 
the probability of occupation for a state in the vacuum is zero, and equation (10) be- 
comes 


v — \ r° 


J = — ) yj M a ij’W* 


4ve 

h / /-«> 


(ID 


The matrix element M ab has been evaluated by Harrison (ref. 4) in the WKB approxi- 
mation, and is given by 


„ 1 - 1/2 - 1/2 

M ah = — P a P h exp 

ab 2 „ a D \ •'x 


■IXA 


(12) 


where x a and x^ are the classical turning points of the motion (see fig. 1), and 




1/2 


(13) 


where e x is the ”x-directed” part of the electronic energy of motion, and V(x) is the 
barrier potential energy. 


V(x) = V D (x) + V t (x) 


(14) 


Calculation of the Matrix Element 

The complete matrix element of the tunneling transition must be calculated not only 
between the initial and final electron states, but between the initial and final vibrational 
states of the molecule as well. If the transition is between the lowest vibrational state 
and a vibrational level of energy hco n , the complete matrix element is <n|M a ^|0) . In 
order to calculate < n | M a b|°>> H is first necessary to evaluate the integral appearing in 
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the exponential of equation (12). Using equations (13) and (14), 



The energy required to excite a molecule from one of its vibrational states to another is 
ordinarily a few tenths of an electron volt (infrared region of the spectrum). The bar- 
rier height V Q (x) - e x , on the other hand, is usually not much smaller than the work 
function of the metal so that V^x) 10 _1 [V o (x) - ej. Setting V Q (x) - e x = <p expand- 
ing the integrand in equation (15), and retaining only the first order term, 



Vj(x) is expected to be important only for the range x < x m , where x m is the point at 
which V Q (x) reaches a maximum (see fig. 1). In this region <p may be treated as a 
constant (i. e. , V Q (x) meets the requirements of the WKB approximation). For x > x 
V i (x) may be ignored completely, but the x dependence of <p m should be considered. 
Hence, the integration is broken into two parts: 



while the second may be written as 
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(19) 


The first integral in equation (19) is a standard expression in the theory of field emis- 
sion from metals and has been evaluated (ref. 6). It is 


-2 



x 


ReF 




where cp is the work function of the cathode and €p its Fermi energy; t(y) and v(y) 
are given below, and have been tabulated by Burgess, Kroemer, and Houston (ref. 7) for 
typical values of electric field and work function. 


v(y) = 




(20a) 


t(y) = v(y) - (*)y 

\3 / dy 


(20b) 


k 2=.. 2l^V ;y=lS (2 

1 + Vl - y 2 * 

where K(k) and E(k) are the complete elliptic integrals of the first and second kinds, 
respectively. The second integral in equation (19) is simply 



( 21 ) 


since the range of integration is that for which cp m is approximately constant. Col- 
lecting results, 
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( 22 ) 


M ab = ^Pa / 2 ' J i 1/2exp 




g , € x ~ e F 
2 2d 


where 


I 





(23a) 


g = 


-4V2m^ 3 

3KeF 



(23b) 


d = 


KeF 



(23c) 


If P x is set equal to zero in equation (22) (i. e. , there is no interaction), the equation 
reduces to the usual expression obtained for field emission from clean surfaces into 


vacuum. 


according to equation (23a), and P x is typically 3x10 


-31 


The quantity £ =* 10 , 

- 19 A 

Coulomb-meter (10” esu-cm) (ref. 2). Hence, the first term in the exponential factor of 


equation (22) is on the order of 10 
is finally written as 


-2 


Using the expansion for exp(-x) for small x, 


<n|M ab |°> 


1 

2u 



^V /2 


1 




This reduces to 


8 



< n Kbl°> ^ V2 Pb 1/2 (— 1 /2 -^ < n l P xl°> ex P ("f + 

V^m/ \ 


€ x ~ g F 
2d 


( 25 ) 


since P x is the only factor in equation (24) which operates on the wave functions de- 
scribing the vibrational states of the molecule. 


Total Energy Distribution 

The part of the current density arising from the electron- molecule interaction is 
given by 





(26) 


where e T = e - e^, and f(e’) is the Fermi function [exp(e'/kT) + 1]~*. In a measurement 
of the total energy distribution only those electrons which satisfy the minimum energy 
condition (see fig. 2) 


e - Ka, n>^col +e F- eV col 


(27) 



Figure 2. - Schematic representation of field emission energy distribution analyzer. Electrons with energy less 
than € will not reach collector. Collector is located far enough from anode that tunneling into collector is 
negligible. 
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will reach the collector. The quantity ho; n is the excitation energy of a vibrational state 
of the molecule, <P co j is the collector work function, V co ^ is the bias voltage between 
cathode and collector, e is the total energy of an electron incident on the barrier from 
inside the metal, and e p is the cathode Fermi energy. The total energy distribution of 
field emitted electrons is defined as 


P T (e) =- — 

L e de 


(28) 


From equation (25) for the "interaction” part of the current, 


pi (e ) .4* 2ss m p *m 






(29) 


where r^ is the "perpendicular" part cf the energy: e = e x + But at the absolute 
zero of temperature in the absence of an interaction, 



Taking the derivative of equation (30) and substituting it in equation (29) results in 


pl ( £ ) = sAie { : 

fra * 




( n 

P x 

°>l 

\X J 

u2 


f(e’) 


(31) 


which, in terms of the applied collector voltage in equation (27) becomes 

2 


, aJ*™ o/r,\ <nP |o>| 

pi T (V col ) = - ^ s 2 Li ) - f („ + ti w - eV ji 

^m dV col \ m / h 2 


(32) 


The total energy distribution at finite temperatures for the noninteracting part of the cur- 
rent density is simply 
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/ dJ 

P T <V col ) = ^ (diT 2 " t^col ' eV col> 


( 33 ) 


coiy 


and Prp(V co i) for the total field emission current density becomes 


P T< V col) = -! 


dJ. 


,dV 


col ) 


-T ,( « , ec.r eV col ) 

e 


, 8, 2 m 2/ r l\ |<"KI 0 >l 


<P 


m 


m/ h" 


f Kol + w n “ eV col> 


From equation (34) 


(34) 


dP r 


d 2 J 


- P^(V^)f dJ ° 


-1 


dJ, 


‘"'col '"COL T C0 ‘V dV -l, 


dV 


col 


1 df(< ?col " e V col> 1 877 2 m 


dV 


col 


<P- 


m \ my 


x 


l<“KI°>| 2 dl <'?’co L + - eV col> 

dV 


h 


col 


(35) 


At low temperatures the derivative of the Fermi function behaves like a Dirac delta func- 
tion, and equation (35) may be rewritten as 


dP r 


d 2 J 




dJ, 


dJ, 


dV col dV 2 ol T CDl \ dV col, 


VdV 


col , 


~ - eV col> * V s 

e 2 <?m 


x £ 


2 / r i 


<n|P x |°)| 


6 <«’col + ®"n “ eV col> 


m/ 


(36) 


Equations (34) and (36) represent the final result for the case where a single molecule 
with a single upper vibrational energy level of energy Kw n interacts with the tunneling 
current. If N foreign molecules per unit area have been deposited on the cathode, the 
average ’’interaction” current density may be found by integrating equation (26) over the 
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surface area of the cathode and multiplying the result by N. Such an averaging process 
assumes that the molecules do not interact with each other to alter their dipole fields. In 
order for the dipole approximation to remain valid, then, the density N must be low 
enough that an electron emerging from the metal sees the field of a single molecule. For 
the idealized case considered here (an infinite plane representing the metal surface), the 
electric field F and work function tp are constants independent of r r Hence, the aver- 
age '’interaction” current density may be written as 




The lower limit of integration over r ± has been set at some point r Q in the vicinity of 
the molecule. (Near the molecule the approximation for the interaction energy fails be- 
cause the transverse components of the dipole moment no longer cancel. ) The integral 
over rj^ is 



Evaluation at the upper limit may be accomplished by noting that 
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Thus 


r 


2 2 2 2 
1 1 + X a + X m + x a x m 


1/2 




lim -< 
r^— °° 


In 


*1 


2 2 2 2 
[1 + _ a + X JP + 


4 

1/2 




+ 1 + 


2 2 
x m +x a 


2r1 


J 


I = 277ln 


4 2/2 2\ 2 2 

r o + r o \ x m + x ay + x m x a 

‘f^n + r o( x m + x a) 
o o \ in a* / 


!/ 2 „ x 2 + x 2 

+ r 2 + _m a 

o o 


2 2 \ 2 2 
x m +x a] + x m x a 


1/2 


= 27rln/3 


(39) 


The exact value of r 0 is unimportant since it enters only in the logarithm. Inserting 
equation (39) back into equation (37) gives 


< J i> = 


327r 3 e 3 Nm 

h\ 

^ m 



(40) 


A molecule ordinarily possesses more than a single vibrational level above its 
ground state. The expression for the ’’interaction” part of the current density given by 
equation (40) must therefore be summed over all such levels, and P T (V col ) and 
dp T ( V col ) / dV col Anally become 


P T^ V col) = - 


dJ 


(dV 


col/ 


\ f ^col " eV col) 


+ N 


327r 3 m 


h> m ^ 


ln(3^ HP x H 2 «?’oo 1 + "n- eV col> 


(41) 
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and 


■ ^T ( V col) =I1Io_\p t (V co1 )I 


dJ 


,-l 


dJ 


dV 


col 


idV 


col/ 


dV 


col . 


,dV 


col/ 


— H<p, 

2 


col “ e ^coP 


+ N-^iBln/3 


,2 

h cp. 


Z l <n N 0> l 


6 ^col + w -- eV — i) 


n 


col 


m 


n 


( 42 ) 


Calculation of the Collector Current Density 

The total average current density is given by 



The average collector current is obtained by integrating over the range of total energies 
for which the electrons are collected. The sum over k^. may be converted to an integral 
as follows: 



h 


Considering only the zero temperature limit for simplicity, 
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< J col> =^~®tp(-g) 


/ 

rr\ 


^col^F eV col 


e -p\ {* / £j\ 

expl — Ide / exp ( ] 

V d / / \ d/ 

«/n 


de , + N 


647r 4 m 2 e 3 


Yi <p 


m 


x lnp\ l^l^xl 0 )! 2 ex P(-g) f 

Z — i nr 


n 


from which 


^col +Rco n +e F- eV col 


/' 

t/ 0 


exp^-^d, l exp(-l% £l 


(45) 


r 


< J col> = ^^exp(-g)<;(<p col - eV col )exp|-_ i-) + d 


T 


(?) 


1 - exp 


*col eV c ol 

d 


cp 


m 



+ N j ffir m 2 e 3 d exp( _ g)ln/3 ^ |<n|P x |o>| 2 ) ( <p CQl + hco n - eV col )exp 


n 


+ d 


i - expp*- + *^l e _Z£2lj 


(46) 


The quantity d ranges in value from a few hundredths of an electron-volt to a few tenths 
of an electron-volt for stable field emission. Since for most metals of interest as field 
emission cathodes e F is on the order of 5 or 6 electron volts, the linear term in each 
pair of braces may be neglected. (The factor exp(-e F /d.) will almost always be less 
than 10" 4 .) 


DISCUSSION 

Figure 3 is a sketch of <J C0 ]> as a function of collector voltage. The dashed 
curve represents the noninteracting part of the collector current density. As indicated 
in the figure, an increase in current occurs whenever the bias voltage is such that elec- 
trons which have lost energy in the electron-molecule interaction can be collected. The 
threshold voltage for each onset of additional current should be exactly equal to 


15 



1.4 



Figure 3. - Collector current as function of collector volt- 
age. Temperature, 0 K. 


ef^cd + Kw n ), where Kw n is the energy of the vibrational level to which the molecule 
makes a transition (assuming the molecule was initially in its ground state). 

Sketches of the total energy distribution and its derivative are shown in figures 4 
and 5. Each increase in <J C0 ]> is accompanied by a corresponding increase in Pi^coP 
for the same values of collector voltage. The threshold voltages are marked by the ap- 
pearance of a "spectral line” in ^T^col^^col' ^ temperatures the lines will 

not be perfectly sharp, as shown in the figure, but will be broadened by the appearance of 
a thermal tail on the distribution function. 

An order of magnitude estimate of the ratio of the "interaction" part of the collector 
current density to the noninteracting part is of interest. As indicated in figure 1, 

O _ O 

x ^ 3x10 centimeter and x o ~ 0, 5*10 centimeter. From equation (39), lnj3 ranges 

m ci _ o 

from approximately 0. 7 for r Q = 0 to about 1.4 for r Q = 3x10"° centimeter, and can be 

set equal to 1 for this estimate. A reasonable value for ^ |(n|P x |0)| 2 is 10“ 61 

2 9 - 38 2 2 ^ 

Coulomb -meter (10 (esu) -cm ) (ref. 2). Again, from figure 1, cp m is ordinarily a 
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Total energy distribution, Pj(V co j) f arbitrary scale 



Figure 4. - Total energy distribution from fiqure 3, indicating presence of electron- 
molecule interactions. Temperature, 0 K. 



Bias voltage between cathode and collector, eV C0 | 


Figure 3. - Derivative of total energy distribution of figure 4, showing exis- 
tence of "spectral lines" corresponding to molecular excitations. Temper- 
ature, 0 K. 
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few electron volts, so the ratio <J co ]> j <J col > is ~ 2x10 N. As indicated earlier, 

the expression obtained for the average current density is valid only when the surface 
density of molecules is low. In computing the matrix element M^, it was mentioned 
that the interaction could be assumed to be negligible for distances greater than x . 
Hence, a value for N of 1 molecule per 20(a) 2 is well within the limits required by the 
approximation made previously to obtain <J col > . The final result is that <J co j) j 
<J co i > 0 ~ 10" 2 . This estimate, although of the same order of magnitude as the calcu- 
lated size of the effect in metal- oxide- metal tunneling, seems surprisingly small in view 
of the dramatic effects observed in field emission microscopes when foreign molecules 
are adsorbed on the cathode surface. The reason for the difference is, of course, the 
approximation used for the interaction potential energy V^(x), and the attendant restric- 
tions to assure its validity. Other types of electron- molecule interactions which have 
been considered to explain the larger changes in current density have been discussed in 
the literature (refs. 6, and 8 to 11). None of them, however, have considered the par- 
ticular effect discussed here. Since the dipole potential is essentially a small perturba- 
tion on the barrier potential, it is not unreasonable to expect that the effect will still be 
present even when other, larger interactions are included. 


CONCLUSION 

The present calculation indicates that electron- molecular dipole interactions should 
give rise to structure in field emission currents, just as they do in tunneling currents in 
metal-oxide-metal sandwiches with impurities trapped in the oxide layer. The magnitude 
of the effect for low surface coverage of foreign molecules is estimated to be of the order 
of 1 percent. The first derivative of the total energy distribution of the emitted electrons 
is directly proportional to the dipole spectral weight function of the molecular species on 
the cathode surface, and should yield a series of "lines” corresponding to each vibra- 
tional transition. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 30, 1969, 

129-02-05-14-22. 
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APPENDIX - ELECTRIC FIELD OF A DIPOLE AT A 
METAL-VACUUM INTERFACE 


From figure 6 the field at a point P in the vicinity of the charge distribution shown 
is 



(47) 



Figure 6. - Dipole and its image at a conducting plane. 


From trigonometry 



+ V 


- 2r^Z cos 6 


and if 


« 



7 




21 


cos 




(48) 


(49) 
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Similarly, 


r\-3 



1 + — cos 9 


(50) 


Since rg = r ^+1, and = r ^ - l ' 


E « 



r l “ ^ r l + 2 1 + — COS 0j 

r l 


i \ 3 


- (rj - Z) 11 cos ff\ 
r l 


(51) 


from which 


E 


3? r i \ Q /- 3Zr’ 

I- l 1 cos 8 1 - — 1 1 ' 1 cos 8' 


( r ^) 3 


1 


(52) 


But 


r' = i(x + h) + jy + kz 


r 1 = i(x - h) + jy + 


kz 


(53a) 

(53b) 


and 


Z' = -iZ x+ jZy^kZz 


l - 


(54a) 

(54b) 


so that near the conducting plane, where 8 ~ 8 ', and for distances such that 
r l ~ r l » h > 


E =* — 


1 L 


2iZ, 


3Zi(2h) 


cos 0 




(55) 
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